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Abstract 

The objective of this paper is to show how the recently proposed method 
by Giusti, Heintz, Morais, Morgenstern, Pardo can be apphed to a 
case of real polynomial equation solving. Our main result concerns the 
problem of finding one representative point for each connected component 
of a real bounded smooth hypersurface. 

The algorithm in yields a method for symbolically solving a zero- 
dimensional polynomial equation system in the afiine (and toric) case. 
Its main feature is the use of adapted data structure: Arithmetical net- 
works and straight-line programs. The algorithm solves any affine zero- 
dimensional equation system in non-uniform sequential time that is poly- 
nomial in the length of the input description and an adequately defined 
affine degree of the equation system. 

Replacing the affine degree of the equation system by a suitably defined 
real degree of certain polar varieties associated to the input equation, 
which describes the hypersurface under consideration, and using straight- 
line program codification of the input and intermediate results, we obtain 
a method for the problem introduced above that is polynomial in the input 
length and the real degree. 
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1 Introduction 

The present article is strongly related to the papers jlO| and Whereas the 
algorithms developed in these references are related to the algebraically closed 
case, here we are concerned with the real case. Finding a real solution of a poly- 
nomial equation f{x) — where / is a polynomial of degree d > 2 with rational 
coefficients in n variables is for practical applications more important than the 
algebraically closed case. Best known complexity bounds for the problem we 
deal with are of the form due to 0, 0, (3|]. Related complexity 

results can be found in 1^, [jl^ . 

Solution methods for the algebraically closed case are not applicable to real 
equation solving normally. The aim of this paper is to show that certain polar 
varieties associated to an affine hypersurface possess a geometric invariant, the 
real degree, which permits an adaptation of the algorithms designed in the pa- 
pers mentioned at the beginning. The algorithms there are of " intrinsic type" , 
which means that they are able to distinguish between the semantical and the 
syntactical character of the input system in order to profit both for the improve- 
ment of the complexity estimates. Both papers jioj and Q show that the affine 
degree of an input system is associated with the complexity when measured in 
terms of the number of arithmetic operations. Whereas the algorithms in ||ic| 
still need algebraic parameters, those proposed in |^ are completely rational. 
We will show that, under smoothness assumptions for the case of finding a 
real zero of a polynomial equation of degree d with rational coefficients and n 
variables, it is possible to design an algorithm of intrinsic type using the same 
data structure, namely straight-line programs without essential divisions and 
rational parameters for codifying the input system, intermediate results and the 
output, and replacing the affine degree by the real degree of the associated polar 
varieties to the input equation. 

The computation model we use will be an arithmetical network (compare to 
pO|). Our main result then consists in the following. There is an arithmetical 
network of size {ndS* L)'-"^^^ with parameters in the field of rational numbers 
which finds a representative real point in every connected component of an affine 
variety given by a non-constant square-free n-variate polynomial f with rational 
coefficients and degree d > 2 {supposing that the affine variety is smooth in 
all real points that are contained in it.). L denotes the size of the straight-line 
program codifying the input and 6* is the real degree associated to f . 
Close complexity results arc the ones following the approach initiated in [^ , 
and further developed in (si), UlII, see also §, §. 



For more details we refer the reader to [hOl and l9(l and the references cited there. 
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2 Polar Varieties and Algorithms 

As usually, let (Q, M and C denote the field of rational, real and complex 
numbers, respectively. The affine n-spaces over these fields are denoted by 
Q", J?" andC", respectively. Further, let (T" be endowed with the Zariski- 
topology, where a closed set consists of all common zeros of a finite number 
of polynomials with coefficients in Q. Let W C C" be a closed subset with 
respect to this topology and let = Ci U ••• U Cs be its decomposition into 
irreducible components with respect to the same topology. Thus W, Ci, . . . ,Cs 
are algebraic subsets ofC". Let 1 < J < s, be arbitrarily fixed and consider the 
irreducible component Cj of W . 

In the following we need the notion of degree of an afiine algebraic variety. 
Let W C C" be an algebraic subset given by a regular sequence /i , • • • fi € 
Q[xi, - ■ ■ , Xn] of degree at most d. If W C (C" ist zero dimensional the degree 
of W, degW, is defined to be the number of points in W (neither multiplicities 
nor points at infinity arc counted). If C C" is of dimension greater than 
zero (i.e. dimW = n — i > 1), then we consider the collection Ai of all affine 
varieties of dimension i given as the solution set inC" of a linear equation system 
Li = 0, • • • , Ln-i = with Lk = Z]"=i akjXj + Uko, au & Q,l < i < n. Let 
A^iv be the subcollection of M formed by all varieties H G M such that the 
affine variety HdW satisfies H nW ^ a.nd dim{H n W) = 0. Then the aflane 
degree of W is defined as max{S\d = deg{H n W), H e Mw}- 

Definition 1 The component Cj is called a real component of W if the real 
variety Cj fl IR" contains a smooth point of Cj . 
If we denote 



then the affine variety W* := (J Cj C C" is called the real part of W. By 

deg*W := degW* = Yl degCj we define the real degree of the set W. 
jei 

Remark 2 Observe that deg*W = holds if and only if the real part W* of 
W is empty. 

Proposition 3 Let f G Q[Xi, ■ ■ ■ , X„] be a non-constant and square-free poly- 
nomial and let V{f) be the set of real zeros of the equaiion f{x) = 0. Assum,e 
V{f) to be bounded. Furthermore, let for every fixed i, < i < n, the real 



I = {j G ]N\1 < j < s, Cj is a real component of W}. 



variety 



Vi := {xGM"\ fix) = 



df{x) 



df{x) 
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be non-empty (and Vq is understood to be V{f)). Suppose the variables to be in 
generic position. Then any point ofVi that is a smooth point ofV(f) is also a 
smooth point of Vi . Moreover, for every such point the Jacobian of the equation 
system f = = ■ ■ ■ = = has maximal rank. 

Proof 

Consider the linear transformation x < — A^^^y, where the new variables are 
y = {Yi, ■ ■ ■ , Yn). Suppose that ^4^'^ is given in the form 

where / and define a unit and a zero matrix, respectively, and 
Oki S R arbitrary if k, I satisfy i -\- 1 < k < n, 1 < I < i. 

The transformation x < — A^^^y defines a linear change of coordinates, since the 

square matrix A^'^^ has full rank. 

In the new coordinates, the variety Vi takes the form 

This transformation defines a map given by 

(^1: ' ' ' ,Yi, ■ ■ ■ , Yn, Oi+i^i, • • • a„i, • • • ai-|-l,i, • • • , an,i) = 

dYi ^ "^■^ay,-'---' dYi ^ 

For the moment let 



Then the Jacobian matrix of is given by 



(i+l)x(n+(n— 
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If aO = (ro, ■ • • ,y„o,aO+i 1, • • -flO^ J belongs to the fibre «'-^(0), where (Fi", • • ■ 
is a point of the hypersurface V{f) and if there is an index j £ {i + 1, • • • , n} 
such that ■^Y' 7^ at this point, then the Jacobian matrix J ($i(a°)) has the 
maximal rank i + 1. ^ 
Suppose now that for all points of V{f) 

dfjy) _ _ dfjy) _ 

and let C := -K" \ = . . . = M£) o}, which is an open set. Then the 

restricted map 

$j : C X iR("-^)* — > 
is transversal to the subvariety {0} in 5?*+^. 

By weak transversality due to Thom/Sard (see e.g. Q) applied to the diagram 

\ i 

one concludes that the set of all A G 5?^"^*^* for which transversality holds is 
dense in K^""*)*. 

Since the hypersurface V{f) is bounded by assumption, there is an open and 
dense set of matrices A such that the corresponding coordinate transformation 
leads to the desired smoothness. □ 

Let / G Q[Xi, ■ ■ ■ ,Xn] be a non-constant squarefree polynomial and let W :— 
{x £(27" I f{x) = 0} be the hypersurface defined by /. Consider the real variety 
V :^W n M"- and suppose: 

• is non-empty and bounded, 

• the gradient of / is different from zero in all points of V 

(i.e. 1^ is a compact smooth hypersurface in _ZR" and / = is its regular 
equation) 

• the variables are in generic position. 

Definition 4 (Polar variety corresponding to a linear space) Let ?, < 

i < n, be arbitrarily fixed. Further, let := {x G C"| Xi^i = ■ • ■ = Xn = 0} 
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be the corresponding linear subspace oJ(U^^. Then, Wi defined to be the Zariski 
closure of 



{xe(F"|/(a;) = 




is called the polar variety of W associated to the linear subspace . 
The corresponding real variety of Wi is denoted by Vi Wi Pi IR^ . 



Remark 5 Because of the hypotheses that F 7^ is a smooth hypersurface and 
that Wi 7^ by the assumptions above, the real variety Vi := Wi PlK", < i < 
n, is not empty and by smoothness of V, it has the description 

{Vq is understood to be V.) 

According to Proposition 3, 14 is smooth if the coordinates are chosen to be in 
generic position. Definition 4 of a polar variety is slightly different from the one 
introduced by Le/Teissier |23]. 



Theorem 6 Let f G Q[Xi, ■ ■ ■ ,Xn] be a non-constant squarefree polynomial 
and let W :— {x gC"| f{x) = 0} be the corresponding hypersurface. Further, 
let V :— W n 5?" be a non-empty, smooth, and bounded hypersurface in M" 
whose regular equation is given by f — 0. Assume the variables Xi, • • • , Xn to 
be generic. Finally, for every i, < i < n, let the polar varieties Wi of W 
corresponding to the subspace X^ be defined as above. Then it holds : 



V C Wo, with Wo = W if and only if f and A ^J^^ 
coprime. 



2 

are 



• Wi is a non-empty equidimensional affine variety of dimension n— + 
that is smooth in all its points that are smooth points of W , 

• the real part W* of the polar variety Wi coincides with the Zariski closure 
in€'^ of 
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for any j, i < j < n the ideal 

df df 



as 



is radical. 



Proof: 

Let i,0 < i < n, he arbitrarily fixed. The first item is obvious since Wq is tlie 
union of all irreducible components of W on which A does not vanish identically. 
Then Wi is non-empty by the assumptions. The sequence /, • • • , of 
polynomials of (Q[Xi, . . . ,Xn] forms a local regular sequence with respect to 
the smooth points of W, since the afline varieties {x eC"|/(a;) = ^qxi = • • • = 
= 0} and {x e(Z7"|^^ = 0} are transversal for any fc,0 < A; < i - 1, by 

the generic choice of the coordinates, and hence the sequence /. • • • . 
yields a local complete intersection with respect to the same points. This implies 
that Wi are equidimensional and dima^Wi = n — (i + 1) holds. We observe that 
every smooth point of W, is a smooth point of W, which completes the proof of 
the second item. 

The Zariski closure of Vi is contained in W* , which is a simple consequence of 

the smoothness of Vi. One obtains the reverse inclusion as follows. Let x* G W* 
be an arbitrary point, and let C^* be an irreducible component of containing 
this point, and Cj* fl Vi ^ 0. Then 

n — i — 1 = dim][i{Cji, fl Vi) = dimjRRiCj* fl Vi) = 

= dim[:R{{Cj» n Vi)') < dirrKcCj^. = n — i — 1, 

where R{-) and ( )' denote the corresponding sets of smooth points contained 
in (•) and the associated complexification, respectively. Therefore, c?iTria7(Cj* fl 
Vi)' = dirridjCj^, = n — i — 1 and, hence, Cj* = {Cj* fl Vi)' , and the latter set is 
contained in the Zariski closure of Vi. 
We define the non-empty afBne algebraic set 



Let j,i < j < n, be arbitrarily fixed. Then one finds a smooth point x* 
in Wi such that ^* ^® fixed in that way. The hypcrsurface 

W = {x e C"|/(.t) = 0} contains smooth point, too. Consider the 

local ring Ow,x' of x* on the hypersurface W. (This is the ring of germs of 
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functions on W that are regular at x* . The local ring Ow,x* is obtained by di- 
viding the ring €[Xi, . . . , Xn] of polynomials by the principal ideal (/), which 
defines W as an affine variety, and then by localizing at the maximal ideal 
{Xi - X'l, . . . ,Xn~ X*), of the point x* = {X^, . .. , X*) considered as a single 
point affine variety.) Using now arguments from Commutative Algebra and Al- 
gebraic Geometry, see e.g. Brodmann j^], one arrives at the fact that Ow,x' is 
an integral regular local ring. 

The integrality of Ow,x* implies that there is a uniquely determined irreducible 
component Y oiW containing the smooth point x* and locally this component 
corresponds to the zero ideal of Ow,x* 7 which is radical. Since the two varieties 
Wi nY and WHY coincide locally, the variety Wi HY corresponds locally to 
the same ideal. Thus, the desired radicality is shown. This completes the proof. 

□ 



Remark 7 If one localizes with respect to the function A(x) — ^ {~B7c~) 
then one obtains, in the same way as shown in the proof above, that the ideal 

dXi''"' dXi''^ 

is also radical. 



Remark 8 Under the assumptions of Theorem 6, for any i, < z < n, 
we observe the following relations between the different non-empty varieties 
introduced up to now. 

V, cv, v^cw* dW.c 

where V is the considered real hypersurface, Vi defined as in Remark 5, Wi the 
polar variety due to Definition 4, W* its real part according to Definition 1, 
and Wi the affine variety introduced in the proof of Theorem 6. With respect 
to Theorem 6 our settings and assumptions imply that n — i ~ \ = dimojWi = 
dimtiWi = dimajW* = dimjpjVi holds. By our smoothness assumption and the 
generic choice of the variables we have for the respective sets of smooth points 
(denoted as before by R{-)) 

V, = RiV,) c RiW,) c RiW,) c R{W), 

where W is the affine hypersurface. 
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For the following we use the notations as before, fix an i arbitrarily, < i < n, 
denote by S* the real degree of the polar variety Wi (compare with Definition 
1, by smoothness one has that the real degree of the polar variety Wi is equal 
to the real degree of the affine variety Wi), put 6* :~ max{(5^,|0 < fc < z} and 
let d :— degf. Finally, we write for shortness r := n — i — 1. 

We say that the variables Xi , . . . , X„ are in Noether position with respect to 
a variety {/i = • • • = = 0} in C", /i, . . . , £ (Q[Xi, . . . , X^], if, for each 
r < k < n, there exists a polynomial of Q[Xi, . . . , Xr, Xk] that is monic in Xk 
and vanishes on {/i = • • • = = 0}. 

Then one can state the next, technical lemma according to where the 

second reference is important in order to ensure that the occurring straight-line 
programs use parameters in (Q only. 

Lemma 9 Let the assumptions of Theorem 6 be satisfied. Further, suppose 
that the polynomials f, , • ■ ■ , ■§x' ^ Q [-^i i ■ ■ • i ^n] are given by a straight- 
line program (3 in (Q [Xi , . . . , Xn] without essential divisions, and let L be the 
size of p. Then there is an arithmetical network with parameters in Q that 
constructs the following items from the input (3 

• a regular matrix of (5"^" given by its elements that transforms the vari- 
ables Xi, . . . , Xn into new ones Yi, . . . ,Yn 

• a non-zero linear form U G (Q[Yr+i, • • . , Yn] 

• a division-free straight-line program 7 in Q[Yi, . . . ,Yr,U] that represents 
non-zero polynomials g G (Q[Yi, . . . , Yr] and q,pi, . . . ,pn £ (Q[Yi, . . . ,Yr,U] 

These items have the following properties: 

(i) The variables Yi, . . . ,Yn are in Noether position with respect to the variety 
W*_^, the variables Yi, . . . , 1^ being free 

(ii) The non-empty open part {W*_j.)g is defined by the ideal {q, gXi —pi, . . . , 
gXn ~ Pn)g in the localization (Q[Xi, . . . , Xn]g . 

(Hi) The polynomial q is monic in u and its degree is equal to 



S*_, = deg* W„-r = dcgVK*„, < S*. 

(iv) max{deg„pfc|l < fc < n} < S*_^, max{degpfc|l < fc < n} = {d6*)°^^\ 
degg= 
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(v) The nonscalar size of the straight-line program 7 is given by {sdS* L)^^^K 

The proof of Lemma 9 can be performed in a similar way as in pO| , |^ for estab- 
lishing the algorithm. For the case handled here, in the i-th step one has to apply 

the algorithm to the localized sequence • ■ ■ , as input. The only 

point we have to take care of is the process of cleaning extranious (^-irreducible 
components. Whereas in the proofs of the algorithms we refer to it suffices to 
clean out components lying in a prefixed hypersurface (e.g. components at infin- 
ity), the cleaning process we need here is more subtile. We have to clean out all 
non-real (^-irreducible components that appear during our algorithmic process. 
The idea of doing this is roughly as follows. Due to the generic position of the 
variables Xi,. . . ,Xn all (^-irreducible components of the variety Wn~r can be 
visualized as (5- irreducible factors of the polynomial q{Xi, . . . ,Xr, U). If we 
specialize generically the variables Xi, . . . , X^. to rational values ?7i, . . . , 77^, then 
by Hubert's Irreducibility Theorem (in the version of |22j) the (^-irreducible 
factors of the multivariate polynomial q{rii, . . . ,r]r,U) correspond to the (Q- 
irreducible factors of the one-variate polynomial q{rji, . . . ,rir,U) G Q[U]. In 
order to explain our idea simpler, we assume that we are able to choose our 
specialization of X\ , . . . , Xr into 771 , . . . , 77^ in such a way^that the hyperplanes 
Xi — rji = 0, . . . , Xr — rjr = cut any real component of Wn-r (This condition 
is open in the strong topology and doesn't represent a fundamental restric- 
tion on the correctness of our algorithm. Moreover, our assumption doesn't 
affect the complexity). Under these assumptions the (5- irreducible factors of 
q{Xi, . . . ,Xr, U), which correspond to the real components of Wn-r, reappear 
as (5- irreducible factors of 9(771, . . . ,r]r,U) which contain a real zero. These Q- 
irreducible factors of 9(7/1, . . . ,r]r,U) can be found by a factorization procedure 
and by a real zero test of standard features of polynomial complexity character. 
Multiplying these factors and applying to the result the lifting-fibre process of 
p^ , [H we find the product q* of the (^-irreducible factors of q{Xi, . . . , Xr, U), 
which correspond to the union of the real components of the variety Wn-r, i-e. 
to the real part of Wn-r- The ideal {q* , gXi — pi, . . . , gXn —Pn)g describes the 
localization of the real part of Wn-r at g. All we have pointed out is executible 
in polynomial time if a factorization of univariate polynomials over (Q in poly- 
nomial time is available and if our geometric assumptions on the choice of the 
specialization is satisfied. 

Theorem 10 Let the notations and assumptions he as in Theorem 6. Suppose 
that the polynomial f is given by a straight-line program (3 without essential 
divisions in (Q[Xi, . . . , Xn], and let L he the nonscalar size of (3. Further, let 
5* := deg*W^i, 5* := max{(5*|0 < i < n} he the corresponding real degrees of 
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the polar varieties in question, and let d :— degf. Then there is an arithmetical 
network of size {ndS*L)'^^^^ with parameters in Q which produces, from the input 
(3, the coefficients of a non-zero linear form u e <Q[Xi, . . . and non-zero 
polynomials q,pi, . . . ,Pn G Q[U] showing the following properties: 

1. For any connected component C ofV there is a point S, G C and an element 
T € IR such that qir) = and ^ = (pi(r), . . . ,Pn (''')) 

2. deg{q) = 6*„_, < 6* 

3. max{deg(pj)|l <i<n} < S*_i. 
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